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Abstract 

The cusps of the caustics of any gravitational lens model can be 
classified into positive and negative ones. This distinction lies on 
the parity of the images involved in the creation/destruction of pairs 
occurring when a source crosses a caustic in a cusp. In this paper, we 
generalize the former definition of the sign of the cusps. Then we apply 
it to the binary lens. We demonstrate that the cusps on the axis joining 
the two lenses are positive while the others are negative. To achieve 
our objective, we combine catastrophe theory, usually employed in the 
derivation of the properties of caustics, with perturbative methods, in 
order to simplify calculations and get readable results. Extensions to 
multiple lenses are also considered. 
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I. Introduction 



The application of catastrophe theory in the study of the lens mapping rep- 
resents a considerable step in the understanding of the critical behaviour 
in gravitational lensing [|IJ. Particularly interesting is the classification of 
singularities through these methods, relying on the evaluation of intrinsic 
mathematical quantities 0. A complete treatment of these topics can be 
found in ||. 

The methods of catastrophe theory have proved very powerful in the study 
of the caustic structure of the binary lens where it has been employed to 
find the transitions between different topologies, previously studied in the 
equal mass case ||. Erdl & Schneider derived the positions of the cusps 
and showed that, in these transitions, pairs of cusps are created or destroyed 
passing through beak-to-beak singularities (4|. Witt & Petters used com- 
plex notation to study the singularities of the binary lens with an additional 
shear field and continuously distributed matter. 

Since the binary lens is one of the most important lens models, a detailed 
study of its singularities can help in the interpretation of the physical be- 
haviour of this system and gain information about features that cannot be 
calculated analytically. For example, the study of the amplification map near 
the cusps, performed in general by Schneider & Weiss |7|], can provide very 
useful information on the amplification of images in some regimes. Several 
results on cusp counting in multiple lens systems can be found in @, ||. 

A particular problem is the creation of images in the neighbourhood of 
cusps. When a source crosses a caustic in a fold singularity, two images 
of opposite parities are created. The creation of these images happens in 
a different way when the crossing occurs at a cusp. In this case, one pre- 
existing image changes parity and two new images of the same parity are 
created. As the parity of the two new images is the same of the first before 
the crossing, the sum of the final parities equals the parity of the single image 
before the crossing. 

The parity of the original image involved in this process is a characteristic 
property of the cusp, called sign [0. Positive and negative cusps clearly 
behave in an opposite way, but also the lens mapping in their neighbourhood 
is influenced in different ways. The sign of a cusp can be determined by a 
detailed study of the analytical form of the caustic, through the evaluation 
of the fundamental quantities of catastrophe theory. 

In this work, we first give an intrinsic definition of the sign of the cusps 



1 



and then we use it to determine the signs of the cusps in the caustics of 
the binary lens. Instead of dealing with the involved exact formulae, we 
prefer to prove our assertions in some particular cases where perturbative 



approximations are available |T(| n| and then extend our results by means 
of continuity arguments. 

In Sect. II, we review the principal steps in the description of a cusp by 
catastrophe theory to state our definition. Sect. Ill contains the body of the 
calculation of the signs of the cusps in binary lenses. Sect. IV contains some 
considerations about the extensions to multiple lenses. 

II. Cusps in catastrophe theory 

As usual, we introduce the Einstein radius of a reference mass M : 
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We indicate the coordinates in the lens plane normalized to R% by x = 
(xx] 22) and the coordinates in the source plane by y = (yi,y2)- All masses 
are measured in terms of M . The matter density normalized to the critical 
density 

S cr = C2D ° S (2) 

is k (x). 

The Fermat potential of a given distribution of matter is: 

(x, y ) = I ( x - y) 2 - - / dV« (x') In |x - x'| . (3) 

2 7T J 

The lens equation is obtained by taking the gradient of this potential: 

V x 0(x,y) = O. (4) 

This equation can be written in the form of an application from the lens 
plane to the source plane: 

y = y(x). (5) 
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Given a source position y, the x's solving the equation (^) are called 
images. 

The local properties of the lens mapping can be studied through its Ja- 
cobian matrix 

where 7i = | (0n — ^22) an d 72 = 0i2 = 4>2i- We use the notation <p i = 
The trace of the Jacobian matrix is 

Tr J = 2 (1 - k) . (7) 

Many properties of the lens mapping can be studied through the deter- 
minant D of the Jacobian matrix 

D = (1-k) 2 -( 1 1 + 1 1). (8) 

The sign of this determinant, evaluated at the position of an image, is called 
parity. Images with negative parity are characterized by a reversed handed- 
ness. In the points where D = 0, the lens mapping is not invertible. These 
points are arranged in smooth curves called critical curves and their images 



through the lens mapping (EH are called caustics [12 



If a parameterization of the critical curve is available in the form: 

x = x c (6) , (9) 
its tangent vector can be written as: 

dx c 

T = ^- < 10 > 

For the caustic, we have the natural parameterization: 

y = y c (0)=y(x c (0)), (n) 

leading to a tangent vector: 

dy c 



d# 



JT = (T- V)y(x)| x=x c (e) . (12) 
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The caustic singularities can be classified into different types, including 
folds, cuspoids, umbilics, beak-to-beaks and lips ||. 

The zero order singularity is the fold, that is a simple critical point, 
satisfying: 

D = 0, (13) 

not belonging to other special classes. 

What we are interested into is the cuspoid sequence. A cusp point is a 
point belonging to the critical curve (thus satisfying condition (0)), such 
that the vector tangent to the caustic vanishes: 



dy c 

w - »■ < 14 > 

Its image on the caustic is the cusp. 

The cuspoid of second order is the swallowtail, coming up when the con- 
dition: 

dV 

^4- = (15) 
d6 2 V ; 



is added to the conditions (O) and (14). The vanishing of higher order 



derivatives leads to singularities of higher order. These definitions do not 
depend on reparameterizations of the critical curve, since for a cuspoid of 
order p we require all derivatives up to the p th order to vanish simultaneously. 

To analyse in detail the lensing behaviour near a cusp, it is convenient 
to choose the coordinates so as to have the cusp point in the origin of the 
lens plane and the cusp in the origin of the source plane. With an opportune 
rotation, the lens mapping can be generally written in the form 

yi = cxi + 2^ x 2 (16a) 

U2 = hx\%2 + axl, (16b) 

where a, b, c are coefficients depending on the derivatives of the Fermat 
potential evaluated at the cusp point |§. 

The determinant and the trace of the Jacobian, in this frame, are 

D = bcx 1 + (3ac - b 2 ) x\ (17) 
Tr J = c. (18) 
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The critical curve is found solving the equation D = 0. We parameterize 
it in this way: 



b 2 - 3ac 

a£ = 0. (19b) 



= ;__^ 2 (19a) 
oc 



The caustic is 



h 2 — ?nr 

yi = 3^P» 2 (20a) 



The sign of the quantity lac — b 2 defines the sign of the cusp . We can 
guess that the behaviour of the caustic changes in an evident way according 
to this sign. To see these differences, one can solve the lens mapping explicitly 
in the neighbouhood of a cusp point. A detailed discussion of this topic can 
be found in Refs. [jj. In Fig. [I], we briefly recall the behaviour of the lens 
mapping near positive (top row) and negative cusps (bottom row). When 
a source crosses the caustic in a positive cusp, a positive image crosses a 
critical curve, thus changing its parity. Moreover, two positive images come 
up from the crossing point in the direction tangential to the critical curve. 
For negative cusps, all parities are reversed. Then a negative image changes 
its parity and a pair of negative images arises. Of course, the sum of the 
parities is unchanged in both processes. 

Now we want to give a definition for the sign of the cusp that can be 
used in all frames. With this definition, we will be able to compute easily 
the signs of the cusps for the caustics of any lens mapping without having to 
look for the standard frame for each cusp. 

The vector tangent to the cusp is 

dy c _ ( 3^0 \ 

~d9~{ 3 b ^6 2 J ■ { } 

Its derivative, evaluated at the cusp is 
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We can notice, recalling (|HS|), that a cusp with a vanishing sign is a swallow- 
tail. 

The vector 

is everywhere orthogonal to the critical curve. In the standard frame, it reads 

N = ( 2(*-L) x J • (24) 
At the cusp point, N is parallel to the second order derivative of the caustic 



N-^|r, (25) 



Then, the quantity 

d 2 y c 
~d9 2 

evaluated at the cusp point, is 3 (2ac — b 2 ). So its sign coincides with the sign 
of the cusp. Moreover, the expression (|25|) is invariant for translations and 
rotations of the coordinate frame and for reparameterizations of the critical 
curve (as long as it is evaluated at the cusp point). Then it is natural to take 
Eq. ( p5|) as the general definition for the sign of the cusp, valid not only in the 
standard frame but for any choice of the coordinates and parameterizations 
of the critical curve. 

In what follows, we will use this definition to determine the signs of the 
cusps in the binary lens. 



III. Signs of the cusps in binary lenses 

The critical curves of binary lenses can be classified according to two param- 
eters: the mass ratio of the two lenses and their separation. We indicate the 
two masses (always measured in terms of the reference mass Mo) by m\ and 
m 2 . The total mass is M = mi + m 2 . For each mass ratio, three possible 
topologies are present depending on the separation a between the lenses. In 
Fig. |2], we show them for a system of two equal masses. 

( J\,/2_ a 4 N ) 3 

For close binaries, i.e. when ^ 27a8 '- > m\mi @, three caustics are 
present: the one at the centre of the system has four cusps, while the other 
two have three cusps (Fig. ||a). 
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Wide binary systems, characterized by the condition a 2 > ym 1 ^ 3 + m^ 3 J , 

have two caustics with four cusps each (Fig. |]c). 

Intermediate binaries have one caustic with six cusps (Fig. |2|b). 

If we change the parameters of the system continuously, we can follow 
each cusp unambiguously, since their positions are continuous functions of 
the parameters. Also the quantity (p5]) , giving the sign of the cusp, is a 
continuous function of the parameters. So, remembering that the two vectors 
in the scalar product of Eq. (p5|) are parallel at the cusp, a change of sign 
in a cusp can only occur through a higher order singularity, according to 
Eq. flI5p. But no higher order cuspoid is present in binary lenses, so each 
cusp bears its sign unaltered as the parameters change, even when transitions 
between different topologies occur (if the cusp is not directly involved in these 
transitions). 

These observations allow us to introduce an overall labelling for the cusps 
as we have done in Fig. |2|. As we are interested in the study of the sign of 
the cusps, we have given the same label to the cusps which can be obtained 
by reflection on one or both axes in the equal-masses case, as they obviously 
have the same sign. For this reason we have two cusps of type a, four of type 
b and so on. 

In total, twelve cusps can be distinguished in the binary lens: two of them 
(indicated by e in Fig. Qc) are only present in wide binaries; four of them 
(indicated by c and d in Fig. only exist in close binaries. The other six are 
always present. The disappearance of cusps c and d and of cusps e occurs 
through beak-to-beak singularities at the critical values of the separation 
between the masses, leading to the intermediate topology 0. 

Since the sign of the cusp is a characteristic not depending on the partic- 
ular choice of the parameters, we can lead our study in the simplest cases and 
then extend our results in general, exploiting the continuity of the quantities 
involved in our calculations. Of course, a natural simplifying choice can be 
the equal-mass case, as concerns the mass ratio. So, from now on we put 
mi — ^2 — M/2. 

Now, we could deal with the exact formulae for the caustics, given in 
Refs. 0, [|, H , but we would rather simplify our treatment further, reducing 
to the cases where perturbative results are employable. In this way, by very 
few steps, we can derive properties of the general binary lens from the study 
of its particular cases. 

In Ref. |13[] approximations are given for the caustics of the binary lens 
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in the form of multipole expansions in its extreme cases: very close binary 
systems (a <C y/~M), very wide binary systems (a ^> y/~M) and planetary 
systems (m2 <C mi). Perturbative expansions to the first significant orders 



for multiple lenses were derived by Bozza JT0[ in the same cases. In the 
following subsections, we use the expressions in refs. |10| , |TT|1 to analyze the 
signs of the cusps in representative cases. Deriving these signs for all types 
of cusps identified in Fig. ||, our demonstration will be concluded. 

A. Cusps in wide binary caustics 

Let's begin from the caustics formed by the two lenses when they are very 
far each other. Choosing the mass of each single lens as the reference mass, 
i.e. M/2 = 1, the starting hypothesis is 

a > 1. (26) 

These caustics are the results of the deformations on the two separated 
Einstein rings of each mass induced by the other. We put the first lens at 
the origin and the second at the position x 2 = (a; 0). Thanks to Eq. (pB|) , 
we shall expand all our objects in series of the perturbative parameter - to 
the first non trivial order. The lens equation becomes: 

x i 1 x i \ 

Vi = xx- —2 H (27a) 

X a Or 

V2 = (27b) 
x a 



The determinant of the Jacobian matrix is 

^ = i- 7 4- 2( l7 i) - ( 28 ) 

|x| cr |x| 

The solution of the equation D = with D given by Eq. (B3) is the 
critical curve of the first mass [[L(| 



cos^(l + ^) 

x c (0) = I I (29) 

sin£(l + ^) 
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The caustic is 



1 , 2cos 3 f 

y c (0) = I I • (so) 

2 sin 3 6 
a 2 

Now we can build the objects involved in the classification of the cuspoids. 
The vector tangent to the caustic 



dy c 3 sin 29 
~d9 ~ aP~ 



cos 9 



sin i 



(31) 



vanishes at 9 = mr /2 for n = 0, 1, 2, 3. The cusp obtained for n = is on 
the Xi-axis, directed towards the other mass. Looking at Fig. |2|c, where the 
caustic of the first mass is the left one, it can be identified as a cusp of type 
e. The cusp for n = 2 is in the opposite direction and is a cusp a. The other 
two are cusps of type b. 

Now, let's calculate the signs of these cusps. The second derivative of the 
caustic is 



^2yC g / cos 9 (— 1 + 3 cos 20) 
sin 9(1 + 3 cos 29) 



d9 



(32) 



and the vector N, defined in Eq. (E3|) is: 



-2 cos 9 + ± (5 cos 9 + cos 39) 
N = | a | (33) 

-2 sin 9 + (-5 cos 9 + cos 36) 



Finally, the quantity (p5|) evaluated at the four cusps is 

12 



N ■ 



c\9 2 



0=7171-/2 



-ir-. (34) 



Then, the cusps a and e are positive, while the cusps b are negative. 
Having established the signs of these three types of cusps, the only cusps to 
be studied are of type c and d, to be found in the analysis of close binary 
systems. 
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B. Cusps in close binary main caustic 



We calculate now the signs of the cusps of the main caustic of close binary 
systems. This is necessary to obtain the sign of the cusps of type c. However, 
we also get a confirmation for cusps of type a. Now we choose the origin to 
be at the center between the two masses and consider the total mass as the 
reference mass: M = 1, x x = [— |; 0), x 2 = (|; 0). We consider the limit 



a < 1. 



(35) 



By virtue of this hypothesis, all objects will be expanded in powers of a. 
Stopping at the first significant order, the lens equation becomes 

a 2 xi (xl — 3^2) 



V2 



Xi 



x 2 - 



|x| 

•i'2 



+ 



4|x| 
a 2 X2 [x\ - 



'ix\ / 



4|x| 



(36a) 
(36b) 



The Jacobian determinant is 



D = 1 



2|x| 



The main critical curve is given by the expression [K| 



V 



sin^ 1 



/ 



The caustic is 



y c (0) 



f cos 3 e 



-\sm s 9 



(37) 



(38) 



(39) 



The calculation proceeds in a way analogous to the previous case. 
The vector tangent to the caustic 



dy^ 
d# 



3a 2 sin 29 



cos 9 



sin 9 



(40) 
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vanishes at 9 = nn/2 for n = 0, 1, 2, 3. The cusps obtained for n = 0, 2 are 
on the xi-axis, thus being of type a. The other two are cusps of type c. 
The second derivative is 



d 2 y c 
~d9 2 



3a 2 



cos 9 + 3 cos 3$ 
— sin 9 + 3 sin 36* 



(41) 



and the vector N is 



N 



-2cos# 



(cos 9 — 3 cos 39) 



-2 sin 9 - ^f- (sin 9 + 3 sin 39) 
Finally, the quantity fl25|) evaluated at the four cusps is: 



d 2 y c 



N- 



d9< 



3a 2 



=nn/2 



(42) 



(43) 



Then, besides the confirmation of the positivity of cusps a, we find that 
the cusps c are negative. 



C. Cusps in close binary secondary caustics 

Finally, we need the sign of the cusps of type d. They arise in close binary 
systems and belong to secondary caustics. Studying the cusps of these caus- 
tics, we also have the confirmation of the negativity of cusps b. We retain 
the same choices of the previous subsection for the parameters of the system. 



The critical curve we investigate is the upper secondary one |TT 



x c (9) = | | (44) 

a I a 3 sin 9 
4 

corresponding to the bottom secondary caustic 



a 3 



(2cos# - sin 29) 

y c (0)=\ | (45) 

_ a + f + Y(2sinfl-cos2fl ) 
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The vector tangent to the caustic 



d9 



a" 
T 



cos 29 + sin 9 



sin 29 + cos 9 



(46) 



vanishes at 9 = 2mr /3 + 7r/2 for n = 0,1, 2. The cusp obtained for n = is 
on the £2-axis, then it is the one of type d. The other two are the cusps of 
type b. 

The second order derivative is: 

cos 9 - 2 sin 29 



d 2 y c 



and the vector N is 



a 
T 



N 



(47) 



- sin 9 + 2 cos 29 



A cos 8 



4 sin 8 



Finally, the quantity (|25|) evaluated at the four cusps is: 

d 2 y c 



N- 



A9 2 



(48) 



(49) 



= -3 Vn. 

0=2n7r/3+7r/2 

All of the cusps are negative, not only the ones of type 6, already studied, 
but also the cusp d. 

In conclusion, in binary lenses, we have that the cusps on xi-axis (a and 
e) are positive and the others are negative. 



IV. Extensions to multiple systems 



In refs. |T0| , [LI], the expressions for the caustics were derived not only for 
binary systems but for a general multiple lens. Then, by further continuity 
arguments, we can extend our results even to some special multiple lenses. 
In a first approximation [ID|, the caustic of a lens very far from the other 



components of the distribution of mass is a Chang-Refsdal caustic 
with shear: 



14, 15 
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\ 



E 

i=2 



rrii sin {2ip i 



+ 



E 

i=2 



rrii cos (2<£>J 



(50) 
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along the direction 



ix) — — arctan 
Y 2 



E 

8=2 



mi sin(2<y3j) 

7i 



E 

,i=2 



mi cos(2i/?j) 
3 



(51) 



In these two expressions, the positions of the masses producing the shear are 
given in polar coordinates p i7 

In Sect. 3.1, we have analyzed the case of a binary system where only the 
mass rri2 produced the shear field on the first mass. This field was directed 
towards m 2 . We found that the cusps along this direction were positive and 
the other two were negative. Then, if we let a third mass m 3 arise from zero 
to some finite value (also placed far from the mass m 1; in order to remain in 
the approximation of wide multiple lenses), the shear would be modified, but 
the shape of the caustic would remain the same. The positions of the cusps 
change continuously and then we can still say that the cusps on the shear 
axis are positive and the other two are negative as long as the topology is 
unaltered. 

Regarding close multiple lenses, we can follow the same procedure. In this 
case, the caustic formed by masses having mutual distances much smaller 
than the total Einstein radius is a quadrupole caustic [ lOfl , with moment 



Q = M 
oriented along 



\ 



i=l 



+ 



mi P 2 i Sil1 ( 2< ^ 



i=l 



(52) 



<p = arctan 



Q-M£m ip 2 cos (2i Pi ) 
i=i 

n 

M E m ipl sin (2<pi) 

i=l 



(53) 



In the binary case, examined in Sect. 3.2, (p = 0. The two cusps on 
this axis were positive and the other two were negative. Again, turning on 
additional masses from zero to some finite value, the orientation changes 
continuously and the signs of the cusps are preserved as long as the topology 
remains the same. Then, the cusps on the orientation axis are positive and 
the other two are negative. 
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For secondary caustics, in close multiple systems, we start referring to 
planetary systems. If we have a star surrounded by one planet placed inside 
the Einstein ring of the star, two secondary critical curves are formed very 
close to the planet [[R|]. According to the arguments in Sect. 3.3, the resulting 
secondary caustics have three negative cusps each. Now, if we add more 
planets inside the Einstein ring of the star, each of them will be accompanied 
by two small secondary critical curves. Since, in a first approximation, each 
couple of secondary critical curves formed by one planet is not affected by 
the presence of the other planets |16|, we deduce that these caustics have 



the same properties of the binary secondary caustics (obtained putting the 
masses of the other planets to zero): they have only negative cusps. Now, 
we let the masses of all the planets increase continuously from their small 
values to their definitive values, in order to obtain a generic multiple system 
with distances smaller than the global Einstein radius. The positions of the 
secondary critical curves will change, but their shape (and the signs of the 
cusps) will remain the same as long as they remain simple critical curves, 
i.e. as long as they do not meet each other to give multiple critical curves 



11 . Yet, even in this case, the transitions can be achieved by suppressing 



two cusps through a beak-to-beak singularity JTTJ] . The surviving cusps will 
anyway preserve their sign. In conclusion, all secondary caustics (both simple 
and multiple) have negative cusps. A direct calculation using the formulae 
given in |l]J confirms this statement derived here by continuity arguments. 



We see that the approach used in this specific problem (calculations in 
perturbative cases and extensions by continuity arguments) can be a powerful 
tool in order to gain information about situations where there is poor analyt- 
ical direct knowledge. In the study of caustics of gravitational lens models, 
this methodology can surely provide a great help for the comprehension of 
the topologies actually realized in physical systems. 
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Figure 1: The lens mapping near cusps. On the right column, we have 
the caustics corresponding to the critical curves on the left. The first row 
represents the situation for a positive cusp. The bottom row shows the 
situation for a negative cusp. The sign of D is indicated in the left figures. 
The circles on the left figures are the images corresponding to the source 
positions in the right. When the source is outside of the caustic (empty 
circles) there is one image, while when it is inside (filled circles) three images 
are formed. Their parity is discussed in the text. 
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Figure 2: Topologies of the caustics of the equal-mass binary lens, (a) Close 
binary, (b) intermediate binary, (c) wide binary. The letters identify the 
different kinds of cusps. 
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